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ABSTRACT 

Using three newly identified galaxy clusters at 2 ~ 1 (photometric redshift) we measure the evolution 
of the galaxies within clusters from high redshift to the present day by studying the growth of the 
red cluster sequence. The clusters are located in the Spitzer Infrared Array Camera (IRAC) Dark 
Field, an extremely deep mid-infrared survey near the north ecliptic pole with photometry in 18 total 
bands from X-ray through far-IR. Two of the candidate clusters are additionally detected as extended 
emission in matching Chandra data in the survey area allowing us to measure their masses to be 
-^^500 — 6.2 ± 1.0 X 10^'^ and 3.6 ± 1.1 x 10^'^ M©. For all three clusters we create a composite color 
magnitude diagram in rest-frame B-K using our deep HST and Spitzer imaging. By comparing the 
fraction of low luminosity member galaxies on the composite red sequence with the corresponding 
population in local clusters at z = 0.1 taken from the COSMOS survey, we examine the effect of a 
galaxy's mass on its evolution. We find a deficit of faint galaxies on the red sequence in our 2; ~ 1 
clusters which implies that more massive galaxies have evolved in clusters faster than less massive 
galaxies, and that the less massive galaxies are still forming stars in clusters such that they have not 
yet settled onto the red sequence. 

Subject headings: galaxies: clusters — galaxies: evolution — galaxies: photometry — cosmology: 
observations 



1. INTRODUCTION 

The redshift range from z = 1 to the present 
day is a particularly dynamic epoch in the history of 
groups and clusters as evidenced by the evolution of the 
morphology-density relation and increa sing fraction of 
blue galaxies with increa- sing redshift (Ca pak et al.ll2007l : 
iButcher fc Oemleii |1984|) . Interestingly, cluster ellipti- 
cals at z ~ 1 already have a narrow dis tribution of red 
color s (the red cluster sequenc e (RCS); iBlakeslee et al.l 
120031 Ivan Dokkum et all [2OOII ). There is some debate 
about the mechanism by which these cluster galaxies ar- 
rive onto the red sequence. It is difficult to distinguish 
whether these red ellipticals all formed their stars and did 
their merging at z > 3, then stopped f orming stars whe n 
they entered the cluster environment (jFord et al.ll2004f ): 
or if they are the product of the merging of gas-poor sys- 
tems which do not produce star formation ()van DokkumI 
[2OO5I) . 

We investigate whether the red population is still in 
the process of forming at z = 1 or if indeed assembly has 
already finished at higher redshift by studying the pres- 
ence of the faint end of the RCS at z = 1 and comparing 
it to the present epoch. We measure the ratio of faint 
to bright RCS galaxies in a sample of three z ~ 1 clus- 
ters from the IRAC Dark Field (described below) . These 
clusters have the benefit of extremely deep 3.6 /xm data 
which allows us to study the faint end of the luminos- 



Electronic address: jkrickOcaltech.edu 

^ Spitzer Science Center, MS 220-6, California Institute of Tech- 
nology, Jet Propulsion Laboratory, Pasadena, CA 91125, USA 

^ Large Binocular Telescope Observatory, University of Arizona, 
Tucson, AZ 85721, USA 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge MA 02138, USA 

* Jet Propulsion Laboratory, California Institute of Technology, 
Pasadena, CA, 91109, USA 



ity function at rest frame near-IR, which traces the peak 
of the spectral energy distribution in galaxies. A con- 
firmed deficit of faint galaxies on the RCS would imply 
that more massive galaxies have evolved in clusters faster 
than less massive galaxies. A constant fraction of faint 
red galaxies between z = 1 and the present would require 
a formation mechanism where galaxies of all masses have 
already joined the red sequence at redshifts higher than 
one. There is evidence that the fain t end of the RCS is 
not completely in p lace by z=0.8 (De Lucia et al.l [20071 : 
iKovama et al.l[2007l and references therein) , although at 
least some clusters at th ese redshifts ap pear to have com- 
plete RCSs to M*-h3.5 (IAndreonll200l . 

These questions are ideally addressed with deep in- 
frared surveys of clusters at high redshift. In the last four 
years deep and wide area surveys in the mid-infrared us- 
ing the Spitzer Space Telescope have substantially opened 
a new window on galaxy and star formation at < z < 3. 
Spitzer now routinely produces imaging of large fields to 
higher resolution and fainter depths than previously pos- 
sible. IRAC, the mid-infrared camera on-bo ard Spitzer, 
takes images at 3.6, 4.5, 5.8, and 8.0 ^m (|Fazio et all 
|2004| ). The shorter wavelengths provide a direct mea- 
surement of the stellar content of galaxies at redshifts as 
high as three. The longer wavelength channels sample 
emission from polycyclic aromatic hydrocarbons (PAHs) 
in low redshift galaxies, as well as direct thermal emis- 
sion from hot dust. The deepest such survey is the 
dark current calibration field for the mid-infrared cam- 
era, commonly known as the "IRAC Dark Field" . These 
deep mid-IR data are supplemented by additional 14 
band photometry including Palomar u',g',r',i', HST/ ACS 
F814W, MMT z', Palomar J, H,& K , Akari 11 & 15 /im, 
Spitzer MIPS 24 & 70 ^m, and Chandra ACIS-I imaging. 
This extensive, multiwavelength dataset gives us the dis- 
tinctly unique opportunity to study the development of 
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the red sequence in galaxy clusters at redshift one. The 
second paper on galaxy clusters in this series will discuss 
the role of star formation in the evolution of clusters at 
z = 1 by examining the Spitzer 24 /im data in conjunc- 
tion with the morphological information from the HST 
ACS dataset. 

This paper is structured in the following manner. In 
Sj2] & fj3] we discuss the data and derived photometric 
redshift determination. Details of the cluster search and 
cluster properties are presented in 2] In iJ5] & [JHl we 
present the color magnitude diagrams of the candidate 
clusters and results of the faint-to bright ratios of red 
sequence galaxies. In Sj7]we discuss the implications for 
the evolution of cluster galaxies. Throughout this paper 
we use Ho = 70km/s/Mpc, f^A/ = 0.3, f^A = 0.7. With 
this cosmology, the luminosity distance at z=l is 6607 
Mpc, but the angular diameter distance is a factor of 
(1 -t- z)"^ less, or 1652 Mpc. All photometry is quoted in 
the AB magnitude system. 

2. OBSERVATIONS & DATA REDUCTION 

2.1. The IRAC Dark Field 

The survey region is the IRAC Dark Field, centered at 
approximately 17h40m +69d. The field is located a few 
degrees from the north ecliptic pole (NEP) in a region 
which is darker than the actual pole and is in the Spitzer 
continuous viewing zone so that it can be observed any 
time IRAC is powered on for observing. These observ- 
ing periods are called instrument "campaigns" , and occur 
roughly once every three to four weeks and last for about 
a week. Sets of long exposure frames are taken on the 
Dark Field at least twice during each campaign totaling 
roughly four hours of integration time per campaign, and 
these data are used to derive dark current/bias frames for 
each channel. The dark frames are used by the pipeline in 
a manner similar to "median sky" calibrations as taken in 
ground-based near-infrared observing to produce the Ba- 
sic Calibrated Data (BCD) for all science observations. 
Each set of dark calibration observations collects roughly 
two hours of integration time at the longest exposure 
times in each channel. 

The resulting observations are unique in several ways. 
The Dark Field lies near the lowest possible region of 
zodiacal background, the primary contributor to the in- 
frared background at these wavelengths, and as such is in 
the region where the greatest sensitivity can be achieved 
in the least amount of time. The area was also cho- 
sen specifically to be free of bright stars and very ex- 
tended galaxies, which allows clean imaging to very great 
depth. The observations are done at many position an- 
gles (which are a function of time of observation) leading 
to a more uniform final psf. Finally, because the cal- 
ibration data are taken directly after anneals, they are 
more free of artifacts than ordinary guest observer (GO) 
observations. Over the course of the mission, the obser- 
vations have filled in roughly uniformly a region 20' in 
diameter. This has created the deepest ever mid-IR sur- 
vey, exceeding the depth of the deepest planned regular 
Spitzer surveys over several times their area. Further- 
more, this is the only field for which a 5-1-year baseline 
of mid-IR periodic observations is expected. 

The IRAC data are complemented by imaging data in 
14 other bands with facilities including Palomar, MMT, 
HST, Akari, Spitzer MIPS, and Chandra ACIS-I. Al- 



though the entire dark field is > 20' in diameter, because 
of spacecraft dynamics the central ~ 15' is significantly 
deeper and freer of artifacts. Therefore, it is this area 
which we have matched with the additional observations. 
The entire dataset will be presented in detail in a future 
paper (Krick et al, in prep). For completeness we briefly 
discuss here the Spitzer IRAC, HST ACS, and Chandra 
ACIS datasets as they are the most critical to this work. 
All space-based datasets are publicly available through 
their respective archives. 

2.2. Spitzer /IRAC 

This work is based on a preliminary combination of 75 
hours of IRAC imaging, which is «30% of the expected 
depth not including a possible warm mission. Even these 
75 hours go well into the confusion limit of IRAC and so 
the additional exposure time will not add a lot of sensi- 
tivity. The Basic Calibrated Data (detector image level) 
product produced by the Spitzer Science Center was fur- 
ther reduced using a modified yersion of th e pipeline de- 
veloped for the SWIRE survey (jSurace et al. 2005 ). This 
pipeline primarily corrects image artifacts and forces the 
images onto a constant background (necessitated by the 
continuously changing zodiacal background as seen from 
Spitzer). The data were coadded onto a regularized 0.6" 
grid using the MOPEX software developed by the Spitzer 
Science Center. 

Experiments with DAOPHOT demonstrate that 
nearly all extragalactic sources are marginally resolved 
by IRAC, particularly at the shorter wavelengths, and 
hence point source fitting is inappropriate. Instead, pho- 
tometry is done using the high spatial resolution ACS 
data as priors for determining the appropriate aperture 
shape for extracting the Spitzer data. We do this by 
first running source detection and photometric extrac- 
tion on the coadded IRAC images using a matched fil- 
ter algorithm with image backgrounds determined using 
the mesh background estimator in SExtractor (Bertin 
et al. 1995) . This catalog is merged with the HST 
ACS catalog. For every object in that catalog ,if the 
object is detected in ACS then we use the ACS shape 
parameters to determine the elliptical aperture size for 
the IRAC images. ACS shape parameters are determined 
by SExtractor on isophotal object profiles after deblend- 
ing, such that each ACS pixel can only be assigned to 
one object (or the background). For objects which are 
not detected in ACS, but which are detected in IRAC, 
we simply use the original IRAC SExtractor photometry. 
Because of the larger IRAC beam, we impose a minimum 
semi-major axis radius of 2". In all cases aperture cor- 
rections are computed individually from PSF's provided 
by the Spitzer Science Center based on the aperture sizes 
and shapes used for photometry. 

Final aperture photometry was performed using cus- 
tom extraction software written in IDL and based on the 
APER and MASK_ELLIPSE routines with the shape in- 
formation from SExtractor, from either ACS or IRAC 
as described above, using local backgrounds. Because 
we use local backgrounds, the measured fluxes of objects 
near the confusion limit should have a larger scatter than 
those non-confused objects, but will on average be the 
correct flux. This will have the effect of adding scatter 
to the CMD described below, but will not cause trends 
of movement for the faint objects in color space. Addi- 
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tionally, this will not effect the photometric redshifts, as 
it will likely shift all IRAC points up or down, but not 
relative to each other. The final 95% completeness limit 
at 3.6 /^m is 0.2 /iJy or 25.7 AB magnitude as calculated 
h'om a number count diagram by measuring where the 
number of observed objects drops below 95% of the ex- 
pected number of objects, where the expected numbers 
are calculated by fitting a straight line to the brighter 
flux number counts. 

2.3. HST/ACS 

The HST observations consist of 50 orbits with the 
ACS comprising 25 separate pointings, all with the 
F814W filter (observed I-band). Within each pointing 
eight dithered images were taken for cosmic ray rejec- 
tion and to cover the gap between the two ACS CCDs. 
The ACS pipeline calacs was used for basic reduction of 
the images. Special attention was paid to bias subtrac- 
tion, image registration, and mosaicing. Pipeline bias 
subtraction was insufficient because it does not measure 
the bias level individually from each of the four ampli- 
fiers used by ACS. We make this correction ourselves 
by subtracting the mean value of the best fit Gaussian 
to the background distribution in each quadrant. Due 
to distortions in the images, registration and mosaicing 
was performed with a combination of IRAF's tweakshifts, 
multidrizzle, and SWarp v. 2. 16.0 from Terapix. The ac- 
tual task of mosaicing the final image was complicated by 
the large image sizes. The single combined mosaic image 
is 1.7GB and reading in all 200 images (160Mb each) for 
combination is impossible for most software packages. 

The final combined ACS image is ^ 15' diameter co- 
incident with the deepest part of the IRAC Dark Field 
and is made with the native 0.05" per pixel resolution. 
Photometry was performed in a standard manner with 
SExtractor. The 3(t detection limit for point sources is 
/ = 28.6(AB). The area common to both IRAC and 
ACS contains ^ 51, 000 detected sources. 

All cluster galaxies are detected in the mid-infrared 
data, and the ACS data are used as priors for extraction 
of the mid-IR photometry. Because the cluster galax- 
ies are detected in the optical ACS images, and cluster 
membership is derived specifically based on the optical 
data, it is the completeness of the optical data thats sets 
the fundamental detection limits. Thus our ability to ex- 
amine the faint end of the cluster (rest frame) K-band 
luminosity function is ultimately limited by the optical 
data, not the infrared data. 

2.4. Chandra/ ACIS-I 

The field-of-view of the Chandra ACIS-I is 17' (for the 
central four chips) , well-matched to the deepest central 
area of the Dark Field. As a result, only a single point- 
ing was required. The lOOks observation was broken into 
three separate observations at different pointing angles 
to reduce the effect of the gap between chips. The data 
was reduced using the newest version of the Chandra In- 
teractive Analysis of Observations software (CIAO 4.0). 
Specific attention was paid to destreaking, bad pixels, 
and background flares. The task merge_all was then used 
to combine event files from the three observations into a 
final event file and image. The final combined Chandra 
image has 0.5" per pixel resolution. Blind pointing is 



expected to be 1". Aperture photometry for extended 
sources is done manually (see §4.1|) . 



3. PHOTOMETRIC REDSHIFTS 

The combined IRAC and ACS catalog contains over 
50,000 objects which makes acquisition of spectroscopic 
redshifts impractical. Even confirmation spectroscopy of 
red galaxies at z = 1 in our three candidate clusters will 
require multiple nights on 8-lOm class telescopes. In lieu 
of spectroscopy we use our extensive multi-wavelength, 
broad-band catalog to build spectral energy distributions 
(SEDs) and derive photometric redshifts. These SEDs 
are fit with template spectra derived from galaxies in 
the Spitzer wide ar ea infrared survey survey (SWIRE; 
lPollettaet"aI]l2007t l. Since the SWIRE templates are 
based on Spitzer observations we find them the best 
choice to use as models for this dataset. Photomet- 
ric redshifts are calculated using Hyperz; a chi-squared 
minimization fittin g program including a correction for 
Gal actic reddening ^Bolzonella et al.ll2000tlCalzetti et al.l 
l200dl . We emphasize that in this work photometric red- 
shifts are used only as a first step to find cluster can- 
didates, and not to determine membership within the 
clusters. 

4. CLUSTER SEARCH 

The mid-IR wavelengths of IRAC are well suited to 
find galaxies at redshift one because the stellar peak of 
the spectral energy distribution has red-shifted into those 
bands. We exploit this fact in a search for clusters. 

The first step in the search is to combine the spatial 
information from the 2-D images with photometric red- 
shifts to visually identify clusters of galaxies. We do 
this by plotting the locations of all galaxies in a certain 
photometric redshift range on both the optical and IR 
images to look for clusterings. We step through red- 
shift space in overlapping z=0.2 bins examining each for 
clusterings. The result of this search is 14 cluster can- 
didates at 0.8 < z < 1.4 with average areal densities of 
10.3 galaxies with similar redshifts per square arcminute. 
Next, for each candidate cluster, we examine the redshift 
distributions of all galaxies within 0.013°of the center 
(~ 375 Kpc at these redshifts). We choose this radius, 
at around one third of the virial radius, as the size at 
which the clustering signal appears to be strongest. The 
candidate cluster redshift distributions are compared to 
the average redshift distributions of 50 regions of the 
same area randomly distributed across the field. After 
this comparison, three excellent candidate clusters re- 
main with peaks in their redshift distributions which are 
greater than 2a above that of the comparison fields (see 
Figure[l]& Table[T]). In addition the three best candidate 
clusters show clear over-densities in both the ACS and 
IRAC images (see Figure [2]). While we cannot guaran- 
tee that this search is complete without extensive spec- 
troscopy, we are confident that we have found the more 
massive clusters at redshift one. 

In a shallow, large-area IRAC survev lEisenhardt et al.l 
(|2008[ ) finds roughly two high-z clusters in an area equiv- 
alent to the IRAC dark field. The DEEP2 survey finds 
seven spectroscopically confirmed groups and clusters 
with 0.75 < z < 1.03 and an upper m ass limit of 
1 X 1 0^'^M(r) in a similar area to our survey (|Fang et al.l 
|2007|) . Our finding of three candidate clusters is in agree- 
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ment with the number of clusters in these other surveys. 
Two of the candidate clusters have centers within 1.8' 
of each other. If clusters were randomly distributed, 
from Monte Carlo simulations we would expect to find a 
close pair like this one 2-50% of the time depending on 
the areal density (here using the IRAC shallow survey 
and DEEP2 respectively). However, from both hierar- 
chical simulations and cluster surveys we do not expect 
clusters to be randomly distributed, instead clusters are 
connected by filament s and are h ighly correlated even 
at redshift one tBrodwin et al.l 1200 71. If our two clusters 
are truly close in all three dimensions then they must lie 
along a filament. Our photometric redshifts do not allow 
for the precision needed to know if they indeed are 3-D 
neighbors. 

4.1. X-ray Properties 

Two of the three candidate clusters show extended 
emission in our lOOks Chandra ACIS-I image. This both 
confirms that they are bound clusters because at these 
detected luminosities they can only be clusters, and al- 
lows us to determine their masses (see Figure [3|). We 
are also able to place a limit on the mass of the third, 
Chandra undetected, cluster candidate. 

X-ray luminosities are obtained by doing aperture pho- 
tometry on the merged ACIS-I image with the CIAO task 
dmextract. Aperture sizes were chosen to be 0.9'(0.4Mpc) 
which fully encompasses all of the excess flux coming 
from the clusters but avoids neighboring point sources 
on the Chandra image. Point sources prevent us from 
using background regions directly adjacent to the clus- 
ters, so instead we use the average of 20 background re- 
gions taken from empty regions all over the frame. Net 
counts are 172 ±31 and 72 ±29 (Poisson errors) for the 
clusters respectively. It is possible that the brighter clus- 
ter (on the right in Figure [3]) is contaminated by a point 
source which is effecting this count rate, however with 
such low number of photons we have not attempted to 
separate this possible point source from the cluster and 
background. 

Taking the count rate from the image, we use the 
Portable Interactive Multi-Mission Simulator (PIMMS v 
3.9d) to derive a flux assuming a thermal bremsstrahlung 
model with 4 x 10^°cm~^ of Galactic NH (measured for 
the location of this field), and a temperature of 5KeV. 
This temperature is chosen randomly and is calculated 
more accurately in the following iterative process. Once 
we have calculated a flux an d luminosity of the clu sters, 
we use a combi n ation of the iMaughan et al.l (|2006( l and 
IVikhlinin erall (|2002l ) L^ - T relation based on Chan- 
dra and XMM data on 22 clusters with z > 0.4 to derive 
a more realistic temperature. We then go back and use 
this new temperature to re-calculate the X-ray flux and 
luminosity. The two clusters with detections have lumi- 
nosities of 3.6 ± 0.6 X lO'*^ and 1.6 ± 0.7 x lO''^ erg/s in 
the 0.5-2 Kev band. 

4.2. Mass 

To derive the mass from the X -ray luminosity f or the 
two detected clusters, we use the iMaughaiil (|2OO70 L^ — 
Afaoo relation based on 34 clusters with 0.5 < 2: < 1.3. 
Derived masses are M500 — 6.2 ± 1.4 x 10^^ and 3.6 ± 
1.4 X 10^'^ Mq. Since the second cluster is a very low 
count detection, we take this mass to be the upper limit 



to the mass of the third, Chandra undetected, cluster. 
Quoted errors only represent error in the measurements, 
including uncertainties in redshift, but do not take into 
account error in the models and should therefore be taken 
as optimistic. 

There is some evidence that using X-ray luminosities 
to estimate mass is unreliable at high redshift in the 
sense that higher redshift clusters have lower x-ray lumi- 
nosit ies for a given mass (jLubin et al.l 120041 : iFang et ahl 
l2007f ). The DEEP2 survey of the aforementioned seven 
clusters at 0.75 < z < 1.0 d etect none of the ir clusters 
in a 200ks Chandra image (jFang et al.ll2007t l. A pos- 
sible reason for this is that clusters are not virialized 
at higher redshifts, an assumption which is necessary to 
use the cluster hot gas to measure mass. Alternatively, 
lAndreon et al.l (2008) have proposed that there are no 
underluminous X-ray clusters and that previous work 
has either not properly measured cluster mass via ve- 
locity dispersion, or have not correctly interpreted their 
X-ray observations. If high-z clusters are underluminous 
in X-rays then there are two relevant implications: 1) our 
high-z clusters have true masses which are higher than 
measured from their Chandra luminosities, and 2) this 
leaves open the possibility that the Chandra undetected 
cluster is a true cluster with mass larger than the limit 
afforded by the current non-detection. The resolution of 
this issue is beyond the scope of this paper. 

As a check of cluster mass we have considered using 
the SDSS relation of optical richn ess to the weak lensin g 
mass of clusters at low redshifts () Johnston et al.ir2007| ). 
This relies on multiple assumptions and relations such 
as that between the number of red sequence galaxies in- 
side of l/i~^Mpc with luminosities greater than 0.4L* 
and r2oo, the radius at which the cluster is 200 times 
the critical density of the universe (i Hansen et al.ll2007l) 
. Furthermore one has to then use the derived value of 
?'200 with the weak- lensing relation to arrive at M2oo- To 
our knowledge none of these relations have been tested 
at high redshifts and we therefore choose not to make 
this calculation. 

The measured mass of these clusters is consis- 
tent with the masses of the, to date, roughly 35 
published groups and clu sters with c onfirmed red- 
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Stanford et al... 2006: .Fang et al...2007. : .Finoguenov et al.l 
2007HEisenhardt et al.ll2008f) . About two thirds of these 



clusters are detected in multi-wavelength surveys such 
as DEEP2, COSMOS, and the IRAC shallow survey and 
have an average mass of 3.8 x IO^^Mq. The rest are de- 
tected in serendipitous X-ray imaging or are targeted for 
their location around radio galaxies. These other clus- 
ters all have higher masses (> li<^14M0) which is con- 
sistent with it being easier to detect high-mass clusters 
in the X-ray and high mass clusters around radio galax- 
ies. Overall we expect from hierarchical formation that 
clusters at 2; = 1 should be less massive than clusters at 
z = 0. 

5. COLOR MAGNITUDE DIAGRAMS 

Most galaxies in clusters are ellipticals that have ap- 
proximately the same red color regardless of magni- 
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tude and therefore form a red cluster sequence (RCS) 
in a color magnitude diagram (CMD), as long as the 
color is chosen to span the 4000 A break. Past stud- 
ies with HST have shown that th e RCS is in place at 
least by redshift one, if n ot before (jBlakeslee et al.|[2003l : 
Ivan Dokkum et al.l[200lh . We choose to make CMDs of 
the dark field clusters with ACS F814W and Spitzer 
3.6 /im. At z ~ 1 this corresponds to rest-frame B- 
K which does span the Balmer break. We generate a 
composite color magnitude diagram for all three clusters 
including all galaxies within 0.017°(^ 500 Kpc) of the 
candidate cluster centers (Figure [3]). The upper axis is 
plotted as absolute K-band magnitude in order to com- 
pare to a low redshift sample. Because 3.6 /im at z = 1 
corresponds to K-band at z ~ 0.1, we do not have to ap- 
ply a K-correction to the data, which avoids a large set 
of uncertainties. We do make a correction for the red- 
shifting of the bandpass and a correction for the luminos- 
ity evolution of galaxie s over time taken f r om th e stel- 
lar evolution models of Bruzua l fc CharlotI ([2003) . The 
color of RCS is consistent with the predicted colors from 
simple stellar populations with a single burst of star for- 
mation above z=3 and a Salpeter initial mass function 
(|Bruzual fc CharlotI [200l . 

Without benefit of spectroscopy we determine mem- 
bership based on location in the CMD. Pinpointing the 
location of the RCS becomes easier when using a compos- 
ite CMD from three clusters. To determine the location 
of the RCS, the bright end of the composite r ed sequence 
([3.6] < 22.5) is fit with a biweight function f Beers et al.l 
|1990| ) where the slope is fixed to -0.12 (to match the low 
redshift slope below). The faint end is left out of the fit 
because of the larger amount of contamination by field 
galaxies. Members are taken to be those galaxies which 
lie within ±0.5 mag of the RCS. We determine the signif- 
icance of the existence of a red sequence in our CMD by 
doing a monte carlo calculation with 100 realizations of 
the CMDs of randomly selected galaxies in the dataset. 
From this calculation we find that the number of galaxies 
on the true RCS is significant compared to random galax- 
ies at the 6.7(7 level. The histogram of member galaxies 
in the composite cluster is shown in Figure [4] 

Although the RCS is mainly composed of cluster galax- 
ies at a particular redshift, there will be some fore- and 
back-ground galaxies that will contaminate the mem- 
bership count. We statistically subtract fore- and back- 
ground galaxies within the measured RCS. The number 
of contaminating galaxies is calculated by counting the 
average number of galaxies from 50 random regions of 
the same size as the clusters which have colors consistent 
with the measured RCS. The average number of contam- 
inating galaxies within a cluster area (0.017°) and their 
standard deviation are 40.3 ± 7.3 and 16.9 ± 8.0 respec- 
tively for the magnitude bins used below ( — 20 < Mk < 
-18 and -22 < Mk < -20). After this subtraction the 
number of member galaxies remaining in each cluster is 
a few to none in the faint bin, and about twenty each in 
the bright bin. 

We create CMDs for a comparison sample of low red- 
shift clusters using the same process on publically avail- 
able B- and K- band data from the C osmic Evolu- 
tion Survey fCOSMOS: ICaDak et al.ll200l . From COS- 
MOS we use four X-ray confirmed clusters with aver- 



age 0.09 < z < 0.12 (jFinoguenov et al.l l2007t Id's 42, 
58, 113, 140). These clusters have an average mass of 
-^500 = 1.1 X 10^'^Mq as determined from the X-ray tem- 
perature. Figure[5]shows both the CMD of the composite 
COSMOS data from four clusters and the K-band distri- 
bution of RCS member galaxies including a statistical 
correction for fore- and back-ground galaxies taken from 
neighboring COSMOS regions. 

This comparison sample is not ideal since these low 
redshift clusters are roughly the same mass as the 
high-z clusters, and we expect clusters to hierarchi- 
cally gain mass over time through the infall of other 
groups and clusters. Unfortunately, to our knowledge, 
there is no sufficiently deep near-IR imaging on local 
~ 10^14}MQclusters. While it is possible that more mas- 
sive clusters will have fewer faint red galaxies than a less 
massive cluster due to merging, we do not expect this to 
be a large effect. Deep, wide-field, near-IR imaging of 
local clusters is required to study this effect in detail. 

The CMD's from the two samples have clear differ- 
ences. Despite our ultradeep data, the distribution of the 
z = 0.1 galaxies extends to fainter magnitudes than that 
for the clusters at z = 1 (thanks to the approximately 
four magnitudes of surface brightness extinction). We 
plot the deeper low redshift data to Mk = —14, whereas 
for the high redshift dark field clusters we can only plot 
a smaller range of magnitudes extending to Mk = —18. 
Everything fainter than that is below our 5cr detection 
threshold, so we are unable to compare that magnitude 
range with the lower redshift clusters. 

6. RESULTS 

We investigate the evolution of the fraction of faint 
galaxies on the red sequence from redshift one to the 
present. To make the comparison between the two dif- 
ferent redshift samples, we choose a rest-frame K-band 
absolute magnitude range which is both brighter than the 
5cr detection threshold for both datasets, and not too far 
towards the bright end of the luminosity function that ex- 
tremely small number statistics would effect the measure- 
ment. We divide the resulting magnitude range into two 
bins which are two magnitudes wide; —20 < Mk < — 18 
(faint bin), and -22 < Mk < -20 (bright bin). 

The quantity we want to measure is the amount of 
faint galaxies on the RCS dX z — 1 compared to the 
present epoch as an evolutionary measure of the build- 
up of the red sequence. This is calculated as the ratio of 
faint RCS galaxies to bright RCS galaxies. This faint-to- 
bright ratio is shown explicitly for both cluster samples 
in Figure [6] and can also be seen by comparing Figures [4] 
& m We find that the relative number of faint galaxies 
on average in clusters at redshift one is less than in the 
average cluster at redshift 0.1 at the 3cr level. 

The major sources of error in this work are the back- 
ground measurement of field galaxies, the area over which 
cluster membership is determined, the possibility that 
the third cluster is a chance projection, and the size of the 
magnitude bins. Recall that completeness in the IRAC 
data is not a limiting factor to this measurement be- 
cause we use ACS source locations and shapes as a prior 
for doing IRAC photometry. To quantify the error orig- 
inating from the field galaxy subtraction, error bars on 
Figure [B] show one standard deviation on the background 
measurement propagated to the ratio of faint to bright 
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number of galaxies. 

To estimate the effect of changing the area over which 
cluster members are counted we re-calculate the faint-to- 
bright ratio using a radius which is 50% larger than the 
original radius (2.25 times the area). The original radius 
was chosen as a compromise between the relative physi- 
cal sizes of the implied area around the low redshift and 
high redshift clusters. For reference the original 500 Kpc 
radius circle is shown on the IRAC images in Figure [21 
The new, larger radius, ratios are shown as stars in Fig- 
ure [6l The effect of changing radius is negligible. 

We examine the possibility that the third cluster is a 
chance projection, not a cluster, and therefore contami- 
nating this measurement. Spectroscopy is the best way 
to confirm that it is a cluster. Without spectroscopy, we 
re-calculate the ratios without the third cluster, using 
just the two X-ray detected clusters. The same trend of 
the bright galaxies far outnumbering the faint galaxies 
on the red sequence at z=l is recovered, albeit with a 
larger noise contribution. Additionally, we test that any 
two clusters are dominating the signal by removing one 
cluster from the measurement in turn. When we recal- 
culate ratios each time the significance of the signal goes 
down. From this we conclude that all three clusters are 
contributing to the signal. 

Lastly, we consider the possibility that the four mag- 
nitude range used is too large in that it goes too close to 
the very uncertain bright end of the luminosity function 
of the z=0.1 clusters, and too far into a regime where 
photometry is confused at the faint end of the luminos- 
ity function in the z=1.0 clusters. To test our results, 
we re-calculate the ratios with 1.5 magnitude wide bins 
centered at Mk = —20. We find, again, the same trend 
of faint galaxies disappearing from the red sequence pop- 
ulation at high redshift. 

There is a clear deficit of faint galaxies on the red se- 
quence at z=l compared to the current epoch when tak- 
ing into account possible sources of error. 

7. DISCUSSION 

We find that there are fewer faint red cluster galaxies 
at high redshift than low redshift in comparison to the 
number of bright red galaxies. There are two possible ex- 
planations for this. There could be overall less faint clus- 
ter galaxies at z=l at all colors than the present epoch, 
or the faint end of the red sequence is not yet in place, 
and instead those galaxies which will fill the faint end at 
z=Q, are still forming stars at higher redshifts. 

The first explanation, that there are overall, at all col- 
ors, fewer faint galaxies at z=l than at z=0, goes di- 
rectly against the theory of hierarchical formation. We 
expect that clusters at higher redshifts will have more 
faint galaxies than today, and that over time the faint 
galaxies will merge into brighter, more massive galaxies. 
This will have the effect of clusters at higher redshifts 
having steeper faint end slopes of the luminosity func- 
tion than today's flatter slopes, or exactly the opposite 
trend of what would be suggested if overall there were 
less faint galaxies at z = 1. While we do not think this is 
a likely solution, it is interes ting that in a si mulation of 
cluster luminosity functions, iKhochfar et al.l (2007) find 
a slight trend that the faint end slope of clusters does 
steepen from redshift 1 to 0, whereas overall they find 
the hierarchical formation trend of fiattening slope from 



z=6 to 0. This is probably due to noise in their calcula- 
tion. It would be nice to know how the overall faint end 
slope of clusters is evolving from z=l to 0, however we do 
not have a large enough sample to attack this problem. 

It is more likely that the deficit of faint red galaxies 
at high z is due to the fact that those faint galaxies 
are still forming stars, and so have bluer colors at those 
redshifts. This implies that a) bright, massive galaxies 
have already shut off their star formation in clusters by 
z=l, and b) faint, less massive galaxies are still forming 
stars in clusters at z=l. We know that clusters exhibit 
red sequences at higher z, so it is not new that massive 
galaxies have already undergone some process, excited by 
the cluster environment, which keeps them from forming 
stars, and lands them on the red sequence. This is fur- 
ther evidence for the popular theory of 'downsizing', in 
which the more massive galaxies evolve first. However, 
all galaxies in clusters do not follow the same evolution- 
ary processes, instead evolution from the blue cloud to 
the red sequence seems to be mass dependent. Whatever 
process stops star formation in clusters appears has not 
yet happened at z=l to the less massive galaxies. Those 
galaxies are still being allowed to form enough stars to 
stay off of the red sequence. The red population of clus- 
ters is not yet fully in place by z=l. 

We do not compare the specific values of faint-to- 
bright ratios of this work with those in the literature 
as different wavelengths, areas, and definitions of faint 
and bright are used. However the trend for a deficit of 
faint gala xies at high re dshifts is in agre ement with th e 
works of [Pe Lucia et ah (2007); Gilbanke jalJ (|2008l ): 
IStott et al.l (12007 ) andJKovama et all ()2007l ). but in con- 
trast to lAndreonI (|2007l ). The cause of the difference is 
unclear, but we note that we use a much redd er, wider 
filter set at high redshifts (rest-frame B-K) than lAndreonI 
(|20b7) (rest-frame u-g) . We also compare with the same 
rest-frame color at low redshift without applying a k- 
correction. 

Future work in this field requires a large enough sample 
of confirmed clusters at high-z with consistent observa- 
tions well below M* to be able to split the sample on clus- 
ter properties. There are hints that the faint-to-bright 
ratios are cluster mass or richness dependent, but the lit- 
erature shows contradictory trends for the se effects likely 
due to the small samples used to date (jGilbank et al.l 
I2008t iDe Lucia et"alll2007[ ). With this evidence for on- 
going star formation in clusters at z=l, in paper two of 
this series we will examine the star formation rates and 
morphologies of the cluster galaxies using deep SPITZER 
MIPS 24 ^m and HST ACS F814W data. Additionally 
as the community continues to build a larger sample of 
high redshift clusters we will be able to study their prop- 
erties, in particular their suitability fo r dark energy num- 
ber count surveys (jWang et al.ll2004l ). 
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TABLE 1 

LUSTER Characteristics 






Cluster 


ra 
J2000 (dcg) 


dec 
J2000 (deg) 




Ngals 

r < bOOKpc 


^peak 


L^ (0.5-2.0 Kev) 
1 X 10«crg/s 


M500 
1 X 1013Mq 


1 
2 
3 


264.68160 
264.89228 
264.83102 


69.04481 
69.06851 
69.09031 




215 
255 

241 


1.0 ±0.1 
1.0 ±0.1 
1.0 ±0.2 


3.6 ±0.6 

1.6 ±0.7 

< 1.6 ±0.7 


6.2 ±1.4 

3.6 ±1.4 

< 3.6 ±1.1 



Redshift peak and one sigma uncertainty are measured from a gaussian fit to tfie redshift distribution. 
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Fig. 1. — Redshift distributions of the three candidate clusters. The soUd histogram shows the redshifts for all galaxies within ~ 375 Kpc 
of the cluster centers. The dotted line shows the average redshift distribution of the field galaxy population measured in regions with the 
same area. These distributions are based on photometric redshifts. 




Fig. 2. — Color composite images of candidate clusters from ACS F814 (blue), MMT z' (green), and IRAC 3.6 /^m (red) data all smoothed 
to the resolution of the IRAC data. The circle shows the 0.017° (500Kpc) radius used to count the number of galaxies on the red sequence. 
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Fig. 3. — Chandra ACIS-I 0.3 - 7kev image of extended emission in the field. The two detected clusters are encircled with a 500Kpc 
radius circle. The third, undetected cluster location is shown with a dashed line circle. The detected cluster on the right has a higher count 
rate and therefore a higher mass. Specifically the extended emission is associated with the center of the cluster and not with the stars in 
the bottom left of Figure [2]The extended emission just below the cluster on the right corresponds to a cluster at Zpiiat=0.25, and is not 
discussed in this paper. 
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Fig. 4. — Left: Color magnitude diagram of all 711 galaxies within 500 Kpc of the centers of the three clusters at z=l from the IRAC 
dark field. Greyscale indicates the density of galaxies within each bin. Lines represent the RCS ±0.5mag. The top scale is absolute rest 
frame K-band magnitude. The dashed line shows the Sit detection limit of the ACS data which does not come close to our fitted RCS. 
The CMD is plotted only to magnitudes well brighter than the measured 95% completeness limit of the IRAC data. Right: Distribution 
of galaxies along the red sequence for the composite cluster at z=l. Dashed lines delineate the faint and bright magnitude bins used in the 
analysis. 
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Fig. 5. — Color magnitude diagram and K-band distribution for the composite COSMOS cluster at z=0.1. Lines have the same meaning 
as Figure |4] For clarity this CMD only includes those galaxies with photometric redshifts in the range of the clusters, however the actual 
measurement is made in the exact same way as for the clusters at 2: = 1 which is to use the RCS minus a statistical background for 
membership information. 
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Fig. 6. — Ratio of faint-to-bright galaxies on the red sequence in both the samples at high and low z. Squares represent background 
subtracted numbers of member galaxies on the red sequence measured inside of a 500 KPC radius. Stars show the same measurement for 
a 750 Kpc radius. Error bars are only shown on the 500 Kpc points and represent error in the background measurement. 



